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ABSTRACT 

We use 14-year baseline images obtained with the Wide Field Planetary Camera 2 on board the Hubble 
Space telescope to derive a proper motion for one of the Milky Way's most distant dwarf spheroidal compan- 
ions, Leo II, relative to an extragalactic background reference frame. Astrometric measurements are performed 
in the effective point spread function (ePSF) formalism using our own developed code. An astrometric refer- 
ence grid is defined using 3,224 stars that are members of Leo II that are brighter than magnitude 25 in the 
F814W band. We identify 17 compact extra-galactic sources, for which we measure a systemic proper motion 
relative to this stellar reference grid. We derive a proper motion [/x Q ,/x,5]=[+104±l 13,-33±151] /ias yr" 1 for 
Leo II in the heliocentric reference frame. Though marginally detected, the proper motion yields constraints on 
the orbit of Leo II. Given a distance of d~230 Kpc and a heliocentric radial velocity v r = +79 km s _I , and after 
subtraction of the solar motion, our measurement indicates a total orbital motion vrj = 266.1±128.7 km s" 1 in 
the Galactocentric reference frame, with a radial component v, c = 21.5±4.3 km s" 1 and tangential component 
v /c =265.2±129.4 km s" 1 . The small radial component indicates that Leo II either has a low-eccentricity orbit, 
or is currently close to perigalacticon or apogalacticon distance. We see evidence for systematic errors in the 
astrometry of the extragalactic sources which, while close to being point sources, are slightly resolved in the 
HST images. We argue that more extensive observations at later epochs will be necessary to better constrain 
the proper motion of Leo II. We provide a detailed catalog of the stellar and extragalactic sources identified in 
the HST data which should provide a solid early-epoch reference for future astrometric measurements. 
Subject headings: galaxies: dwarf, individual (Leo II); Local Group; proper motions; Galaxy: halo 



1. INTRODUCTION 

The Leo II dwarf spheroidal (dSph) galaxy was origi- 
nally discovered in the cou rse of the Palomar Sky Survey 
dHarrington & Wilson! 119501) . Its great distance from the 
Milky Way (MW) marked it as bei ng sign i ficant, so it w as the 
target of early studies (e.g. lHodgdll962l) . iVogt et al.l (11995b 
found evidence f or dark matter fr o m its 11 kms" 1 velocity 
dispersion, while Mighe ll & Rich! (119961) obtained the first 
HST-based color magnitude diagram (CMD) from WFPC2 
photometry, deriving an age of 9 ± 1 Gyr for its intermediate 
age population. The most recent distance estimates, measured 
from t he tip of the red gia nt branch, place Leo II at 233 ± 15 
Kpc dBellazzini et al.ll2005l) from the MW, and it remains as 
one of the most distant known MW satellites that are candi- 
dates to be bound to th e Galaxy. A recent a nalysis of 17 1 stel- 
lar radial velocities by iKoch et all d2007al) finds no evidence 
that the galaxy is experiencing the effects of tidal disruption, 
and indicates a mass to light ratio of 25-50 (M/L) Q . The 
Leo II dwarf galaxy is of great interest not only as a distant 

1 Based on observations made with the NASA/ESA Hubble Space Tele- 
scope, obtained from the Data Archive at the Space Telescope Science In- 
stitute, which is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS 5-26555. These observa- 
tions are associated with program #1 1 182. 
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3 Department of Physics and Astronomy, University of Leicester, Uni- 
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5 UCLA, Department of Physics and Astronomy, 430 Portola Plaza, Los 
Angeles, CA 90095, USA 

6 Leiden Observatory, Leiden University, PO Box 9513, 2300RA Lei- 
den, the Netherlands 



companio n that offers a const raint on the mass of the Milky 
Way (e.g. Watk ins et al.ll2010T) . but also as one of the nearest 
galaxies that might have conceivably undergone evolution in 
isolation from the effect of the Milky Way. Furthermore, its 
orbital path provides an important feature of the evolutionary 
history of the Leo II dwarf. 

Proper motions of other nearby dSphs have been detected 
and estimated using images from the Hubble Space Tele- 
scope (HST). A proper m otion was determined for the For- 
nax dS ph (d » 140 kpc) bylDinescu et all (120041) . iPiatek et al J 
( 12001). and IPiatek et al.1 (|2007|), for Ursa Minor (d-66 kpc) 
by IPiatek et al.l J2005I), Sculpt or (d~79 kpc ) by IPiatek et all 
(2006), Carina (d~ 100 kpc) IPiatek et all d2002l), and the 
Small Magellanic Cloud (dc±58 kpc) by IKallivaval il et alj 
(2006). In all cases, the orbital integration based on current 
Galactic mass models yields bound orbits, with apogalactic 
distances not exceeding ^150 kpc, although including the 
proper motion uncertainties, one cannot rule out apogalactic 
distances as high as 300 kpc f or Sculptor at the 95% confi- 
dence level (Piatek et al. 2006). However, recent proper mo- 
tion measurements of the Large and Small Magellanic Clouds 
yield space velocities larger than expected, which suggests 
that the two objects may be only marginally bound to the 
Galaxy ( IKallivavalil et al.ll2009l) . 

Due to their relatively large distances, the proper motions 
of dSphs are at the sub-pixel level over the baselines typi- 
cally available from HST archival images (< 20 years). For 
example, at a distance of 100 kpc, a 100 km s" 1 velocity trans- 
lates into a proper motion of only 0.211 mas yr" 1 , which 
yields a net motion of only 0.03 pixels on the WFPC2 cam- 
era over a period of 15 years. Due to the undersampling of 
the point spread function in HST images, the detection of 
such a small astrometric motion requires a special astromet- 
ric reduction procedure. The effective Point Spread Function 
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(ePSF) method, pi oneered by I Anderson & King! ( 120001) and 
IPiatek et al] d2002l) . is well suited to the task. 

We have developed our own reduction software, based on 
the ePSF method, in order to detect and calculate the proper 
motion of Leo II. This paper presents initial results from 
this analysis, which allows us for the first time to investigate 
whether this remote dSph is an actual, bound satellite to the 
MW and, if so, to study in detail its orbital characteristics. 
The archival HST datasets used in our study are described in 
§2. The astrometric reduction is summarized in §3. The iden- 
tification of extragalactic sources, critical in establishing an 
astrometric reference frame, is described in §4. The proper 
motion of Leo II is calculated in §5. The inferred space mo- 
tion and orbit of Leo II is discussed in §6. 

2. HST OBSERVATIONS AND ARCHIVAL DATA 

The dSph Leo II was first observed with the Wide Field and 
Planetary Camera 2 (WFPC2) on May 15, 1994. In addition 
to a pair of shallow 80s exposures in the F555W band, a to- 
tal of eight deep exposures (600 seconds each) were obtained 
in each of the F555W and F814W bands. Dithering patterns 
were not generally used at the time, and the 1994 frames have 
no significant offsets between them. 

Leo II was reobserved on March 19, 2004. This time, only 
eight exposures in the F814W band were obtained, 500-700 
seconds each. The 2004 fields are aligned with the 1994 fields 
to within w 3", and with the same field orientation. The 2004 
exposures followed a 4-point dithering pattern. Frames 1 and 
2 constitute a pair with no offset between each other, as are 
the pairs consisting of frames 3-4, 5-6, and 7-8. The 4-point 
dithering pattern occurs between each of the pairs. 

Leo II was observed again on March 25, 2008, yielding 
eight exposures of 1,100 seconds in each of the F555W and 
F814W bands. This time an 8-point dithering pattern was 
used, so no two images are perfectly aligned, and display 
small offsets between each other. The 2008 frames have the 
same orientation as the 1994 and 2004 frames, but are offset 
by 20" in the direction of Right Ascension. This means that 
on one side of each WF camera frame from 1994/2004, there 
is a band » 200 pixels wide which is not covered by the 2008 
exposures. This leaves an area approximately 4 arcmin 2 , or 
75% of the field of view of the Wide Field Camera chips, 
which was imaged at all three epochs. The center of the 
WFPC2 field of view was offset only w 1' from the center 
of the dwarf spheroidal. W ith Leo II having a core radius of 
2.64' (IColeman et al.ll2007h . the HST observations thus cover 
~ 22% of the area within the core radius of Leo II. 

3. ASTROMETRIC ANALYSIS 

The positions of all sources on the frames were re- 
calculated in the effective point spread function (ePSF) for- 
malism, using a method analogous to the one introduced 
and developed b y Anderson & King (2000) (hereafter AK) 
and IPiatek et"afl (120021) . The method consists of fitting the 
pixel profile of a source [Pk]ij = Pk(i,j) with a local rep- 
resentation of the response of the camera to a point source 
VP* = Xie,j— yk), where (xk,yk) denote the hypothetical 
location of the point source in the (z,/) pixel grid. Instead 
of being modeled by a mathematical function, the functional 
form of vf, the ePSF, is extracted from the data using the mul- 
tiple samplings provided by the large number of point sources 
in the field, because each of the k sources provides samplings 
of ^> at the locations (i-XkJ-yk)- Because the hypothetical 
locations of the point sources (xk,yk) are not known a priori, 
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FIG. 1. — Statistical distribution of the second moment of the distribution 
A (or "sharpness parameter") for sources detected in the Leo II field by the 
WFPC2 camera. A value of A m indicates that the source has a spread 
similar to the local ePSF, which means it is unresolved by the camera. The 
distribution shows a tail of objets with A > 1.4; these are sources which are 
effectively resolved by the camera. 

the functional form of ^ must be determined using an iter- 
ative procedure, where the (Xk,yk) are re-calculated and the 
functional form of \fr re-evaluated a fter each iteration, as de- 
scribed in (And erson & K ing 2000). We have developed our 
own reduction software, which determines the ePSF directly 
out of the data, using an iteration procedure, exactly as pre- 
scribed in AK. Our ePSF is determined separately for each 
of the three epochs, and independently for each band. An 
ePSF is also calculated separately for each of the Wide Field 
Camera frames (WF2, WF3, WF4), and the ePSF is allowed 
to vary continuously across the frame, following its separate 
determination in each of nine sectors across the chip. Since 
our general procedure is in most poin ts identical to the one 
described in Anderson & King! (120001) . we do not repeat the 
details here. 

Once a satisfactory ePSF is determined, an algorithm is 
used to converge the centroid (xk , y^) of each source. The algo- 
rithm differs from a least-square fit in that it uses a weighting 
of all the pixels proportial to the local first derivative of the 
ePSF profile, which gives more weight to the profile edges 
than to the central pixel which often yield little positional in- 
formation. The flux ft of the source is re-evaluated after each 
iteration from an ePSF fit of the stellar profile calculated for 
the current (Xk,yk) estimate. Convergence is attained after 5- 
10 iterations. A x 2 test is then used to evaluate the good- 
ness of fit; sources which differ significantly from the ePSF - 
such as camera artifacts, cosmic ray events, and very extended 
sources (e.g. galaxies) - are rejected. The first moments of the 
residuals Tf, are calculated as: 

T k = Y}.Pk\i-fkm-x k ,j-y k \ (1) 

hi 

and should be «0 assuming is a close estimate of the flux. 
Second moments of the distribution A* are then calculated, 
which offer a quantification of how much intrinsic "spread" 
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TABLE 1 

Stellar members of Leo II defining the astrometric grid 3 



a 


5 








F555W 


F814W 


F555W-F814W 


A 


camera 


(J2000) 


(J2000) 


(mas yr~') 


(mas yr~' ) 


(mag) 


(mag) 


(mag) 




WF... 


loo.353iyo2±U,UUUUUu7 


■ n 1 /10C7CO I r\ (\c\r\c\(\c\i 

+22.1485 /5o±U.UUu00U/ 


n 1 A 1 1 QC 

-0.34±0.36 


u.y/±u.3 7 


25.90 


25.26 


0.63 


0.17 


2 


i o icct/Ito I a a a a aa aa 

loo.355243o±u.(X)u0uuy 


i 11 1 AOf\OCA 1 A A A A AA AA 

+22. 14oUo54±U.0Uu0u0y 


0.51 ±0.36 


a 1 1 a /in 
U.12±0.4y 


25.98 


25.41 


0.56 


-0.29 


2 


168.3553641 ±0.0000001 


i ii 1 c i it 1 a I a aaaaaai 

+22. 1 5 1 27 1 0±0. 0000002 


0.07±0.06 


A A C 1 A 1 A 

0.05±0.10 


22.84 


21.93 


0.91 


0.09 


2 


1 *io iccTont I a aaaaaa^ 

loo.3553oyo±U.0UUU0Uo 


■ 11 1 /(01C1A 1 A AAAAAA^ 

+22. 14y252U±U.UUU00Uo 


A 1 1 I A 1£! 

U.71±0.2o 


A 1C\ 1 A 1 1 

-0. 79±0. 33 


26.02 


25.59 


0.43 


0.01 


2 


i /"o iccoaaa 1 a aaaaaai 

168. 3553999±0. 0000002 


j ii i ci n/"oi 1 a aaaaaa 

+22. 1 5 1 9683±0. 0000004 


0.47 ±0.14 


A AT 1 A 11 

0.07±0.22 


25.17 


24.71 


0.46 


0.14 


2 


1 o ^ c c a f\/~ r\ I a aaaaaao 

168. 3554060±0. 0000003 


• ii 1 /)n 1 A/ 1 1 a aaaaaai 

+22. 1473 1 96±0. 0000002 


A 1 "1 1 A 1 A 

0.11±0.10 


A A 1 I A 1 A 

0.01±0.10 


24.52 


23.73 


0.80 


0.05 


2 


168. 3554 /69±0. 0000006 


+22. 1465200±0. 0000004 


A A1 I A 1 1 

-0.03±(J.21 


0.34±0.21 


25.38 


24.84 


0.55 


0.20 


2 


168. 35554 /8±0. 000000/ 


+22. 1454925±0. 0000005 


A A 1 I A 1 A 

0.0l±0.29 


0.52±0.24 


25.39 


24.77 


0.63 


0.08 


2 


168.355562 / ±0.0000005 


+22. 1502803±0. 0000006 


-0.4l±0. 22 


0.26±0.2y 


25.45 


24.83 


0.62 


0.02 


2 


168.3555798±0.0000009 


+22.1513628±0.0000005 


l.03±0.33 


0.69±0.28 


25.80 


25.20 


0.60 


0.03 


2 


168.3556039±0.0000003 


+22. 1465930±0.0000()02 


0.l4±0.l0 


-0.12±0.11 


23.88 


23.05 


0.82 


0.19 


2 


168.3556093±0.0000008 


+22.1532404±0.0000004 


0.30±0.22 


0.08±0.17 


25.84 


25.30 


0.54 


-0.13 


2 


168.3556196±0.0000009 


+22. 1484238±0.0000()04 


0.48±0.29 


-0.13±0.23 


25.08 


24.55 


0.53 


0.16 


2 


168.3556291±0.0000003 


+22.1536420±0.0000()05 


0.0l±0.l5 


-1.01±0.24 


25.61 


25.06 


0.56 


0.10 


2 


168.3557186±0.0000010 


+22. 14855 11 ±0.0000008 


-0.29±0.44 


-0.75±0.38 


25.79 


25.15 


0.64 


0.26 


2 


168.3557224±0.0000005 


+22. 1 542072±0.0000()04 


-0.2l±0.23 


-0.40±0.22 


25.79 


25.32 


0.47 


0.20 


2 


168.3557244±0.0000002 


+22.1532839±0.0000006 


0.4l±0.l3 


-0.41±0.30 


25.26 


24.75 


0.51 


0.09 


2 


168.3557257±0.0000005 


+22.1510173±0.0000005 


-0.34±0.20 


0.59±0.25 


25.83 


25.15 


0.68 


-0.07 


2 


168.3557325±0.0000003 


+22.1573483±0.0000005 


0.67±0.l4 


0.08±0.28 


25.96 


25.49 


0.47 


0.18 


2 


168.3557424±0.0000049 


+22. 1484209±0.0000()09 


0.18±0.82 


-0.24±0.62 


26.31 


25.93 


0.38 


0.14 


2 



a This table is available in its entirety in the electronic version of this paper published in the Astrophysical Journal website. The first 20 lines of the table are shown here as a guide to 
the table format. 
Relative proper motions in the defined Leo II stellar grid. 



there is in the source compared to the local ePSF: 

u 

(2) 

The second moment is expected to be ~ only if the source 
has approximately the same spread as the local ePSF — i.e. it 
is unresolved by the camera. But if the source has a profile dif- 
ferent from the ePSF, then could be significantly different 
from 0. This value A* thus works like a "sharpness parame- 
ter", distinguishing between resolved and unresolved sources. 
We find that the vast majority of the stars in the Leo II field 
have -0.4 < A/. < 0.4, which determines the range over which 
sources can be considered point-like; we use those sources 
to define the local astrometric grid. Sources with a second 
moment A^ > 0.4, have a significantly larger spread than ex- 
pected from the local ePSF, and are effectively resolved by the 
camera. This is clear in Fig. 1 where we see a tail of sources 
with larger A^ values. Visual examination of sources with 
large measured values of A^ indeed reveals the sources to be 
noticeably extended compared with most other sources in the 
field. Sources with a very large sharpness factor (A > 3) are 
so broadened compare to the ePSF that they generally fail the 
X 2 test and are rejected by the code; hence the sources which 
are detected by the code but flagged as "resolved" are still 
relatively compact objects. 

The "unresolved" sources are used to build a relative astro- 
metric frame. Astrometric corrections are performed in mul- 
tiple steps, in order to bring the positions measured in all ex- 
posures/epochs into a single common reference grid. Each 
unresolved source is used in defining the grid with a weight 
proportional to its estimated astrometric uncertainty; in ef- 
fect this yields a larger weight on the brighter sources. Cor- 
rections are performed first on individual frames of a given 
epoch/band; these correct for small offsets, rotations, or dila- 
tions between successive images, notably the offsets from the 
dithering patterns and also changes in the field scale due to 
thermal dilation in HST and the WFPC2. Average positions 



are then recalculated for each epoch/band using the mean po- 
sition of all the measurements from that particular set. Addi- 
tional corrections are then performed to bring the calculated 
mean positions from all seven epochs/bands into a common 
master grid. Finally, proper motions are calculated for each 
individual source, by a linear regression of the mean position 
of the source as a function of epoch. 

Astrometric and photometric data for all of the unresolved 
(point-like) sources is provided in Table 1 . The table lists the 
measured positions and relative proper motions of the 3,224 
sources, along with F555W and F814W magnitudes. We also 
list the second moments A of the source profiles. The WF 
camera frame (2-4) in which the source was detected is also 
listed. These stars define our astrometric reference grid. 

4. IDENTIFICATION OF EXTRA-GALACTIC REFERENCE SOURCES 

A CMD of all of the sources in the combined Leo II field 
is shown in Fig. 2. Distinct plots are presented for the re- 
solved (bottom panel) and unresolved (top panel) sources, as 
defined in Fig. 1 . The CMD for the unresolved sources clearly 
displays the asymptotic giant branch, horizontal branch, 
red giant sequ ence, and main sequence expected for the 
dSph (see also iMighell & Rich! [19961 iColeman et afll2007t 
iKomivama et al.l 120071) . A distinct population of blue strag- 
glers is also detected. The CMD of the resolved sources 
tells another story. The giant sequences are indistinguishable. 
However, a significant number of objects are detected on or 
near the main-sequence locus. Because their colors and mag- 
nitudes are consistent with main sequence stars, we conclude 
that these "resolved" sources are most likely visual binaries, 
probably due to crowding in the dense Leo II field. 

The CMD of unresolved sources also displays a distinct 
population of very red objects (F555W-F814W>1.4 mag). 
Many of these are moderately bright and within the range 
(V < 23.5) from which we can expect accurate astrometry. 
Since no comparable sources are detected in the unresolved 
population, we conclude that these are, in fact, extended ob- 
jects, with point spread functions that are locally wider that 
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TABLE 2 

Extra-galactic reference sources in the Leo-II field 



ID 




a 




S 


Ma" 


M<5 a 


F555W 


F814W 


F555W-F814W 


A 


camera 






(J2000) 




(J2000) 


(mas yr~') 


(mas yr~' ) 


(mag) 


(mag) 


(mag) 




WF... 


EGR-01 


168 


,3577273±0. 0000024 


+22. 


1 75 1 675 ±0.0000007 


0.36±0.82 


-0.29±0.37 


25.33 


22.23 


3.10 


0.64 


3 


ECjR-02 


168 


3578068±0. 00000 1 8 


+22. 


1766783±0. 0000010 


-0.03±0.70 


f\ AO 1 A C 1 

0.48±0.51 


26.78 


23.92 


2.85 


0.73 


3 


EGR-03 


168 


.3625 /U5±U.UUUUUU7 


+22. 


1525 70U±U.UUUUUU4 


0.52±0.24 


-U.0o±0.24 


25.47 


23.15 


2.32 


0.53 


2 


EGR-04 


168 


3629439±0.0U0U015 


+22. 


1628222±0. 00000 1 1 


-0.90±0.66 


-0.32±0.60 


25.41 


23.58 


1.84 


1.03 


2 


EGR-05 


168 


.363596 / ±0.0000008 


+22. 


1521 349±0. 0000004 


0.53±0.26 


f\ Ac 1 n ,1 

-0.65±0.24 


25.91 


23.49 


2.42 


0.51 


2 


EGR-06 


1 AO 

loo 


.3O5034 / ±0.000001 1 


+22. 


1 C/^flOl 1 _l_fl flAflAflAfi 

1 5D053 1 ±0.0000009 


0. /o±0.49 


i noj_A /in 
1.0o±0.4y 


OC 

25.24 


n oa 
22.96 


2.2o 


0.73 


2 


EGR-07 


168 


3658434±0.0000004 


+22. 


1526715±0.0000007 


1.09±0.32 


0.37±0.39 


24.27 


21.65 


2.62 


0.59 


2 


EGR-08 


168 


3667410±0.0()00007 


+22. 


1786564±().()00()0()6 


-0.25±0.39 


0.73±0.28 


24.76 


22.33 


2.43 


0.46 


3 


EGR-09 


168 


3665803±0.0000012 


+22. 


1513390±0.0000010 


-1.13±0.50 


0.14±0.52 


25.87 


23.18 


2.69 


1.10 


2 


EGR-10 


168 


3684006±0.0000003 


+22. 


1693588±0.0000008 


-0.46±0.26 


-0.10±0.51 


25.02 


22.58 


2.44 


0.75 


3 


EGR-11 


168 


3718307±0.0000008 


+22. 


1709446±0.0000005 


-0.49±0.36 


-0.62±0.30 


23.82 


21.59 


2.24 


0.42 


3 


EGR-12 


168 


3815968±0.0000001 


+22. 


171 8691 ±0.0000009 


-0.38±0.12 


-0.87±0.47 


25.07 


22.70 


2.37 


0.48 


4 


EGR-13 


168 


3853748±0.0000007 


+22. 


1681213±0.0000006 


0.71±0.34 


0.68±0.34 


24.61 


22.61 


2.01 


0.46 


4 


EGR-14 


168 


3858761±0.0000006 


+22. 


1704062±0.0000014 


-1.23±0.47 


-1.53±0.75 


24.62 


22.57 


2.04 


1.05 


4 


EGR-15 


168 


.3900591 ±0.0000008 


+22. 


1725934±0.0000010 


1.47±0.46 


0.70±0.55 


25.06 


23.09 


1.97 


0.91 


4 


EGR-16 


168 


3910732±0.0000003 


+22. 


1715817±0.0000004 


-1.08±0.16 


0.77±0.29 


23.90 


21.48 


2.42 


0.56 


4 


EGR-17 


168 


.3927718±0.0000006 


+22. 


1666944±0.()00()0()6 


-0.18±0.24 


-0.42±0.42 


24.86 


22.53 


2.32 


0.47 


4 



a Relative proper motions in the defined Leo II stellar grid. 

the fiducial ePSF. This, combined with the fact that they are 
found significantly off the main stellar locus, suggests that 
they are background galaxies. A close examination of the 
sources confirms this impression (Fig. 3). 

We use the CMD to separate probable Leo II members from 
extragalactic background sources. The idea is that most extra- 
galactic objects should fall off the stellar locus. Distant galax- 
ies in particular should be redder than stellar members as a re- 
sult of their intrinsically red populations and redshift. To help 
in the identification of an extragalactic reference set, we par- 
ticularly consider the distribution of resolved sources, which 
are likely to include most of the background galaxies (see 
above). The distribution of resolved source in Fig. 2 is unam- 
biguous. A significant number of resolved objects define a lo- 
cus well to the red of the fiducial stellar locus of Leo II mem- 
bers. Sources with F555W-F814W>1.3 and F555W<25.0 
are indeed very unlikely to be stellar members of Leo II. We 
formally identify as extra-galactic sources all those which fall 
inside the empirical box plotted in Fig. 2, except for one faint 
source with a very large astrometric uncertaintry. 

The resolved (Leo II) sources, which define the local astro- 
metric grid, and the extra-galactic sources, which are used as 
background astrometric reference, are all plotted in Fig. 4 as a 
function of R.A. and declination. Extragalactic reference ob- 
jects are distributed over all three independent camera frames. 
The extragalactic sources are also relatively spread out and 
sample both the centers and edges of each field. 

Astrometric and photometric data for the extra-galactic ref- 
erence objects is provided in Table 2. The table lists the mea- 
sured positions and proper motions of the 17 sources, along 
with F555W and F814W magnitudes. We also list the sec- 
ond moments A of the source profiles (see §3 above), which 
measures how the profiles compare to a point source. The WF 
camera frame (2-4) in which the source was detected is also 
listed. 

5. RELATIVE PROPER MOTION OF LEO II MEMBERS 

The identification of a set of extra-galactic sources makes it 
possible, in principle, to estimate a proper motion for Leo II 
by measuring the systemic proper motion of these extra- 
galactic sources relative to the astrometric reference frame 
defined by the stellar members of Leo II. The absolute proper 



motion of Leo II with respect to the extragalactic rest frame is 
then simply the vector opposite to this relative proper motion 
of extragalactic sources. As this represents the proper motion 
measured from the heliocentric rest frame, the component of 
the solar motion around the Galaxy must the be taken into ac- 
count to determine an absolute proper motion relative to a ref- 
erence frame at rest relative to the Galactic barycenter, from 
which the orbit of Leo II can be calculated. 

First we examine the proper motion in the direction of the X 
and Y pixel positions {pLx, Hy)- These are different from R.A. 
and Dec, because each of the three WF2, WF3, and WF4 
cameras has its X-Y axis oriented differently. The X and Y 
positional accuracy is thus more sensitive to instrumental ef- 
fects (systematic errors in particular), and better sets the astro- 
metric accuracy of our measurements. The proper motion dis- 
tribution is plotted in Fig. 5 (left panels) as a function of mag- 
nitude. Sources brighter that F814W=24.0 — a total of 634 
objects — have mean values (/Ix,iIy) = (-0.23,-0.08) mil- 
lipixels yr" 1 with a dispersion {0^,0-^) = (2.54,2.60) mil- 
lipixels yr" 1 . The low mean values are consistent with these 
brighter sources being part of the astrometric reference grid 
(of which they are a subset). The dispersion values yield 
an estimate of the intrinsic uncertainty in the proper motion 
measurement for individual point sources, which is thus on 
the order of ±0.0025 pixels yr" 1 along both lines (X) and 
columns (Y) on the WFPC2 frames. The fact that the scat- 
ter values are very nearly identical in X and Y suggests that 
charge transfer efficiency (CTE) effects are negligible in this 
context, otherwise one would expect a larger scatter in the 
proper motion along the columns (Y). Sources fainter that 
F814W=24.0 — 2,593 sources — have mean proper motion 
values (px,iIy) = (+0.20, +0.01) millipixels yr" 1 with a disper- 
sion ((T^.a^y) = (6.30, 6. 10) millipixels yr" 1 . The astrometric 
errors are thus a factor 2.5 larger than for the brighter sources; 
the scatter is also observed to increase as the magnitudes get 
fainter. 

To further examine possible CTE effects, we plot fi x a a 
function of X and [i Y as a function of Y (Fig. 5, right panels). 
CTE effects yield systematic offsets between the centroid of 
bright and faint sources at large Y. Degradation of CTE over 
time would increase this offset in the later epochs, resulting in 
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FIG. 2. — Color magnitude diagram of the Leo II field. All stars identified 
and measured with our astrometric code are shown. The upper panel shows 
the stars which remain unresolved after the ePSF fit. The lower panel shows 
all of the resolved objects (sources with a sharpness factor A >1.3). All of the 
unresolved sources produce a CMD consistent with Leo II membership. The 
resolved sources fall along two distinct loci: they are either coincident with 
the low-mass stellar locus of Leo II and hence are likely visual binaries, or 
they have very red colors, which then strongly suggest that they are extended 
extra-galactic sources (galaxies/quasars). Objects enclosed in the dashed-line 
box are formally identified as extragalactic sources. One faint source within 
the extragalatic cut (black dot within the dashed area) was rejected because 
of a large astrometric uncertainty. 

an apparent proper motion fly of the faint sources relative to 
bright ones. In Fig. 5, we plot fix,fiy as individual points for 
the bright sources (F814W<24.0), and plot the mean proper 
motions fix, fly for the faint sources, in 7 bins along the X 
and Y positions (broken lines). For the faint sources, the dis- 
persion <7ji Z ,<7),j, is noted for each bin with errobars. We ob- 
serve that the bright sources do not show any systemic offset 



in their proper motion as a function of X or Y, except for per- 
haps a weak increase in fiy for Y > 600. The faint sources 
have their mean proper motions fix, fly within 3a of 0.0 for 
all X and Y; there is, however, a weak trend of fix increasing 
with X, and fly increasing with Y, which could possibly be 
due to CTE effects. However, the extra-galactic objects are 
also bright sources, and we expect them to behave more like 
the bright stars, which show negligible CTE effects, if any. 

The 17 extra-galactic (resolved) sources have mean motion 
values (fix, fly) = (-0.18,-1.37) millipixels yr" 1 with a disper- 
sion (jjpx = (6.71,8.04) millipixels yr" 1 . Their scatter 
is thus significantly larger than that of the stellar (resolved) 
sources. Additionally the scatter appears to be larger in Y, 
though this may be due to small number statistics. In any case, 
the scatter in the motions is clearly larger for the extragalac- 
tic sources by a factor of about 3. Part of this scatter could 
be due to a net relative proper motion of the extra-galactic 
sources: this is because the WF2, WF3, and WF4 frames have 
their XY grids rotated by 90 degrees; a net offset in (fi a ,fis) 
would result in a scatter in (/ix,/iy), since the sources come 
from different frames. Indeed, the opposite should be true for 
systematic offsets in (/ix,/irX which would conspire to yield 
a large scatter in (fi a ,Hs)', see below. Beside the possibility 
of a net proper motion, the larger scatter in the extra-galactic 
sources could simply be due to their resolved nature. The 
ePSF would not be the best model for those sources, which 
would result in larger astrometric errors. Individual ePSF pro- 
files, tailor-made for eac h source, would yield more a ccurate 
results as suggested by Mahmud & Anderson (2008); how- 
ever such models require large numbers of sampling (> 100) 
of each object individually, whereas the ePSF of point sources 
is built from the sample of all of the point sources on the 
image, requiring fewer exposure frames. Our current dataset 
does not contain nearly enough frames of Leo-II to build such 
individual ePSF profiles. 

Reported in the equatorial frame, the stellar sources in 
Leo II with F814W<24 mag have mean proper motion 
(< fi a >,< fig >)=(-2.7,-7.8) /ias yr" 1 with a dispersion 
{afx a ,ufii)=(21\.9,291 ' .2) /ias yr" 1 . The dispersion yields an 
estimate of the accuracy in the proper motion of invidual 
sources, while the near-zero value of the mean proper motion 
confirms the stability of our astrometric reference frame. The 
distribution of proper motions for these unresolved sources is 
shown in Fig. 6 (dots). 

The proper motion distribution of our 17 selected extra- 
galactic sources is also plotted in Fig. 6, where we also plot the 
estimated internal uncertainties in the proper motion estimates 
for each individual source (errorbars). The unweighted mean 
value is (/T Q ,,/I ( 5)=(+31.9,+25.9) /ias yr" 1 and the points have 
a dispersion (cr/iQ,,<7/i,5)=(790.5,683.9) /ias yr" 1 . Since one 
pixel on the WF cameras is pa lOOmas, a 700/ias yr" 1 motion 
is equivalent to 7 millipixels yr" 1 . Hence the scatter in (fJ, a ,fJ,s) 
is comparable to the scatter in (/ix,/iy). However, the scatter 
in the data points is significantly larger than the estimated sta- 
tistical uncertainties of the individual points, which are ±450 
/ias yr" 1 on average. This indicates that the proper motion 
of the extra-galactic sources suffers from systematic errors, 
which are probably due to the sources being resolved and the 
local ePSF a marginal fit to their profile. The dispersion in the 
data points suggests that the uncertainties on individual mea- 
surements are underestimate by a factor « 1.7. We account 
for these systematic errors by increasing the uncertainties of 
the individual proper motions by the same factor. We then es- 
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FIG. 3. — Finder charts for the extra-galactic sources used as astrometric reference objects. All fields are 14.5" on the side with North up and East left. Most 
sources are associated with an extended diffuse emission which identifies them as background galaxies. 
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FIG. 4. — Position of the astrometric reference stars from Leo II (black dots) 
and reference extra-galactic sources (red triangles). Extra-galactic sources 
are found in each of the WFPC2 wide field camera frames. Note the rectan- 
gular shape of the fields which is due to an offset in the 2008 epoch images. 
The surveyed area covers 0.0011 square degrees near the center of Leo II, 
representing «i22% of the area within the galaxy's core radius. 

timate the systemic proper motion of the extragalactic sources 
by calculating their weighted average. We find the systemic 
proper motion to be (< fi a >,<fis >)=(— 104±1 13,+33±151) 
/ias yr" 1 , where the quoted errors are the 1-er uncertainties. 
The 1-er, 2-er, and 3-cr limits are plotted as ellipses in Fig. 6. 

The proper motion of Leo II in the extra-galactic frame 
is the opposite of this vector. We therefore estimate 
that Leo II has an absolute proper motion (< fi a >,< 
fig >)=(+104±113 ,-33±151) /ias yr" 1 , which is the proper 
motion as observed from the heliocentric rest frame. If Leo II 
had no net motion in the plane of the sky, we would expect 
the dwarf spheroidal to have a net proper motion due to the 
reflex motion of the Sun in the Galactic rest frame; based on 
the distance to Leo II, this reflex motion is predicted to yield 
a relative proper motion (fi a ,fi5) re f[ ex =(+80,+l&6) fias yr -1 . 
This value is denoted by the crosshairs in Fig. 6, and must 
be accounted for in calculating the orbital motion of Leo II. 
Adding up this component due to the reflex motion of the 
Sun yields for Leo II a proper motion (/j ct o,/t<so)=(+184± 1 13, 
+153 ± 151) /ias yr" 1 in the Galactic rest frame. This indi- 
cates that the transverse motion of Leo II is detected at the 
2-er confidence level. 

6. SPACE VELOCITY AND THE ORBIT OF LEO II 

To evaluate the space motion of Leo II and, in particular, 
to constrain the parameters of its orbit, it is useful to report 
the position and motion of Leo II in a reference frame cen- 
tered on, and at rest with respect to, the Galactic center. The 
most straightforward method is to start with the angular posi- 
tion (I, b) and proper motion (/!/ , fib) in the galactic coordinate 
system, and include both the measured distance from the Sun 
(r = 233 ± 15 kpc) and heliocentric radial velocity (v r = +79 
km s" 1 ) of Leo II. The position (r,l,b) and motion (v r ,fii,fib) 



in this spherical, right-handed coordinate system have a Sun- 
centered, Cartesian equivalent (X,Y,Z) and (U,V,W) defined 
as: 

X = r cos/ cos b (3) 

Y = r sin I cos b (4) 
Z=r sinfo (5) 
I/ = — 4.74 r (fi/sinl + fij, cos/sin/? ) + v r cos/cos/? (6) 

V = 4.74 r (/iicosl -fit sinZsinb ) + v r sin/cosb (7) 
W= 4.74 r fi h cosb +v r sinb (8) 

with X, Y, and Z expressed in kpc, and U, V, W expressed in 
km s" 1 . 

Here we recalculate the angular position and distance of 
Leo II in a galactic coordinate system defined by a Cartesian 
set {Xq,Yq,Zq) and (t/ G ,V G ,W G ) such that the system has its 
origin at the Galactic center and, by definition, is also at rest 
with respect to the Galactic center. Positions and motions can 
be transformed to this system from the Sun-centered positions 
and motions following: 



X G =X-7.94 
Y C = Y 
Z C = Z 

U G = U + 10.00 

V G = V + 5.25 + 225.0 

W G = W + 7.17 



(9) 
(10) 

(11) 
(12) 

(13) 
(14) 



where we have a dopted the galactocentric d istance of 7.94 kpc 
for the Sun from lGroenewegen et alJ (120081) . a relative motion 
of the Sun to the local stan dard of rest (U Q ,V0,W )= (+ 10.00, 
+5.25, +7.17) km s" 1 dDehnen & Binnevl Il998ah . and a 
Galactic rotation velocity of 225 km s" 1 at the solar vicin- 
ity. In this system, we calculate for Leo II a position 
(X c ,y G ,Z G )sa(-76.8, -58.2, 214.8) kpc. We convert the 
motion into this Galactic reference frame, using a Monte 
Carlo simulation to evalutate the range of uncertainties, and 
find (C/ g ,Vg,W g )=(100.8±126.6,216.0±156.9, 118.1±49.3) 
km s , for a total Galactocentric space velocity of vgrf = 
266.lil28.7kms" 1 . 

This Galactocentric Cartesian system also has an associated 
spherical coordinate system (r G) / G ,£> G ) and (v rG ,/X£ G ,/i£ G ), 
which are related to the Cartesian system (X G ,F G ,Z G ) and 
(Uc,Vc,Wc) with the same relationships as in Eqs. 3— 
8. In this system, we calculate that Leo II is located at 
(r G ,Z G , Z> G )«(235.5 kpc, 217.132°, +65.839°). Its Galacto- 
centric radial velocity and proper motion is (v ro ,fii G ,fib G ) - 
(21.5±4.3 km s~\ -99.8±147.8^as yr" 1 , 215.6±1 12.8/ias 
yr" 1 ). This proper motion yields a total transverse velocity 
v, G =265.2±129.4 km s" 1 , a large value compared to the es- 
timated radial motion v rc =21.5±4.3 km s" 1 . This yields a 
mostly tangential space velocity, from which we infer that 
Leo II either has a low-eccentricity orbit, or is currently close 
to perigalacticon or apogalacticon distance. 

With these kinematic components at hand we integrated its 
orbit for a variety of Galactic potentials, where the integra- 
tion was carried out 14 Gyr backwards and forward in time. 
All potentials use d for our test cases have disk and bulge 
components (e.g. 



iDehnen & B innev 



Paczvnski 1990: lAllen & Santillanl [19911: 
1998b), which, however, do not contribute 



significantly to the potential given the large distance (in par- 
ticular its large height above the disk) of the dSph. More 
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FIG. 5. — Proper motions fix and fir of the Leo II field sources along the camera X and Y coordinates, for all three WF cameras. Unresolved stellar sources are 
plotted as black dots, while the extra-galactic reference objects as plotted as red triangles (with errorbars). Panels on the left show proper motions as a function 
of the source magnitude in the F814W band. The stability of the reference frame allows for proper motion measurements with a mean accuracy ±4 millipixels 
per year for unresolved sources brighter than F814W=24, but resolved extragalactic sources show a noticably larger scatter. Panels on the right show fix and fiy 
as a function of column number X and row number Y, respectively. Individual points are shown for bright sources (F814W<24.0, both stellar and extra-galactic) 
while the mean values are shown for the fainter sources (broken line) with 1-cr uncertainties noted by errobars. A weak trend at large Y values indicates possible 
effects from charge transfer efficiency. 



importantly, this large distance renders remote satellites like 
Leo II impo rtant tracers of the Galactic halo and its total mass 
distribution (Wilkinson & Evans 1999). This is already obvi- 
ous from plotting the escape velocity of such models at the 
position of Leo II (Fig. 7). 

In this figure, we sho w the best-fit scale-free potential 
from lWatkins et all {2010) with its total mass of (2.7 ± 0.5) x 
10 12 M© that has been established using tracer satellites out to 
larger distances of ^300 kpc, but we note that the same argu- 
ments hold for the variou s (dark) MW halo models found in 
the li t eratur e to date (e.g. Battaglia et al. 2005t iDehnen et ail 
I2005ll2006h . Even if our proper motion measurements were 
significantly in error and assuming that only the radial ve- 
locity component contributed reliably, Leo IPs space velocity 
of VcRF,rad = 266 ±128 km s -1 is still high and similar to the 
value we obtain from using the full kinematic information. 



Thus, it is the large distance and radial velocity which govern 
its derived dynamics. The comparison in Fig 7 implies that 
Leo II is a true, bound satellite to the MW: the local escape ve- 
locity exceeds the space velocity of Leo II by l.Ocr. We note, 
however, that the concept of "escape velocity" for the present 
arguments should be taken with caution, since its definition 
commonly neglects all mass outside of the considered radius. 
Moreover, recent evidence suggests that the rotation velocity 
of the Milky Way may be larger than used in the present work 
and all current G alactic potential models, at 0o ~ 250 km s" 1 
dReid et al.l 120091) . We will address these issues in a future 
work dealing with the implications of Leo EE's boundedness 
for current MW mass models (Koch et al. in prep.). 

In Fig. 8 we show the orbital solution based on the afore- 
mentioned integrations and our fiducial proper motion mea- 
surements. While currently formally "bound" to the MW po- 
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FIG. 6. — Relative proper motion of the sources in the Leo II field. Extra- 
galatic sources are plotted in red and as filled triangles, with errorbars show- 
ing the estimated 1-sigma measurement errors. The concentric circles show 
the lcr, 2cr, and 3cr confidence limits for the systemic proper motion of the 
extragalctic sources in the Leo II reference frame. The crosshairs are centered 
on the expected motion of the extragalactic sources if Leo II were a static tar- 
get (i.e., the expected proper motion caused by the reflex motion of the Sun 
in the extragalactic rest frame). 
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FIG. 7. — Escape velocity as a function of Galactocentric distance for the 
best-fit halo model from Watkins et al. (2010), which adopts a total mass out 
to 300 kpc of 2.7 X 10 i2 Mq. Both under exclusion (triangle) or inclusion 
(cross) of our proper motion measurement it appears likely that Leo II is 
bound to the MW at present. Also shown are the 1, 2, and 3a- contours). 

tential, it is evident that Leo II has come a long way over its re- 
cent past. Currently at pericenter, its "orbital" period implied 
from the calculations is 50 Gyr and, rather than an apocenter, 
it reaches its largest distance from the MW at the limits of our 
integrations, at ^1.8 Mpc. Considering these time scales and 
large distances from the actual Galactic potential, any timing 
argument and orbital solution becomes necessarily unreliable 
and one would need to account for the entire dynamical his- 
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FIG. 8. — Orbital solution from our best velocity measurements in Galactic 
coordinates. Red and blue curves refer to backwards and forward integration, 
respectively. Its current location (blue solid circle) and the position of the 
MW (cross) are also indicated, as well as its position along the orbit at ±(2, 
4, 6, 8, 10) Gyr (small points). 

tory of the Local Group d Peebles et al. ll989HBvrd eTaIlll994t 
iPasetto & Chiosill2007l) . 

In order to better understand Leo IPs orbital characteris- 
tics, we recomputed its parameters in a Monte Carlo sense by 
varying position, velocity, and proper motion within their un- 
certainties. As a result, we find periods shorter than 6.5 (9.1, 
11.6) Gyr in only 5% (10%, 15%) of the realizations. Like- 
wise, in only 4% (6%, 12%) of the cases can we reproduce 
pericenters within 100 (150, 200) kpc, thus bringing Leo II 
into th e realm of the majority of the M W's closer, bound satel- 
lites dGrebel etalJ l2003: Kochl l2009l) . 

Therefore, we conclude that this dSph is rather an isolated 
Local Group satellite that is falling into the MW regions and 
passing its (dark) halo for the first time (e.g., Chapman et al. 
2007; Majewski et al. 2007). In fact, interactions or common 
origins of dynamical interlopers li ke Leo II with o ther Local 
Group systems appears plausible (ISales et al.11 2007). 

7. CONCLUSIONS 

Our astrometric analysis provides the first-time detection 
of a proper motion for the dwarf satellite galaxy Leo II. The 
transverse motion of Leo II is detected, but only after account- 
ing for the solar reflex motion. 

The question of its origin and whether or not Leo II is a 
bound satellite has important implications for our understand- 
ing of the infall and merging of cosmological subhalos and 
the dynamics and structure of the Local Group. On the other 
hand, Watkins et al. (2010; their Fig. 5) estimate that the cu- 
mulative contribution of Leo II to the mass budget of the MW 
is typically less than 8% (thus < ±0.2 x 10 I2 M Q ). Removing 
it from the family of bound MW satellite tracers would thus 
only have a minor impact. 

We can compare the idea that Leo II spent most of its life in 
seclusion from th e MW to the properties of other isolated Lo- 
cal Group dSohs. iMidiell & RW(fl9^ . lDolphid(l2002l) and 
iKoch et al] (|2007b) find evidence of star formation as recent 
as 2 Gyr ago, which is comparable to another, remote MW 
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dSph, Leo I dGallartet al.lll999l) . Tucana (at 900 kpc) and Ce- 
tus (775 kpc from the MW) are associated with neither M3 1 
nor the MW and neither of them shows any eviden ce of past or 
present interactions with th ese massive galaxies dLewis et alJ 
l2007t IFraternali et al1 l2009l) . Despite the lack of any obvious, 
dynamic gas removing agents at those large distances, none 
of these dSphs contains any gas; in particular Tucana appears 
to have had no activ e star formation over the past 8-10 Gyr 
dSaviane et al.ll9 96). This may imply that galaxies like Leo II 
had experienced efficient galactic winds, which is also consis- 
tent with its flat age-metallicit y relation during t he first seven 
or so Gyr after the Big Bang (Koc h et al.1 
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